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Abstract 

Background: Recently, we reported an information density tlieory and an analysis of three- para meter plus shorter 
scan than conventional method (3P+) for the amyloid-binding ligand [^^QPittsburgh compound B (PIB) as an 
example of a non-highly reversible positron emission tomography (PET) ligand. This article describes an extension 
of 3P + analysis to noninvasive '3P++' analysis (3P + plus use of a reference tissue for input function). 

Methods: In 3P++ analysis for [^^QPIB, the cerebellum was used as a reference tissue (negligible specific binding). 
Fifteen healthy subjects (NC) and fifteen Alzheimer's disease (AD) patients participated. The (index of receptor 
density) values were estimated with 40-min PET data and three-parameter reference tissue model and were 
compared with that in 40-min 3P + analysis as well as standard 90-min four-parameter (4P) analysis with arterial 
input function. Simulation studies were performed to explain kj, biases observed in 3P++ analysis. 

Results: Good model fits of 40-min PET data were observed in both reference and target regions-of-interest (ROIs). 
High linear intra-subject (inter-15 ROI) correlations of kj, between 3P+-I- (/-axis) and 3P + (X-axis) analyses were 
shown in one NC (/^ = 0.972 and slope = 0.845) and in one AD {i^ = 0.982, slope = 0.655), whereas inter-subject k^ 
correlations in a target region (left lateral temporal cortex) from 30 subjects (15 NC + 15 AD) were somewhat lower 
(f^ = 0.739 and slope = 0.461). Similar results were shown between 3P-I-+ and 4P analyses: = 0.953 for intra-subject 
k^ in NC, = 0.907 for that in AD and = 0.71 1 for inter-30 subject k^. Simulation studies showed that such lower 
inter-subject k^ correlations and significant negative kj biases were not due to unstableness of 3P-i~i- analysis but 
rather to inter-subject variation of both /c2 (index of brain-to-blood transport) and kj, (not completely negligible) in 
the reference region. 

Conclusions: In [^^CjPIB, the applicability of 3P++ analysis may be restricted to intra-subject comparison such as 
follow-up studies. The 3P-i~i- method itself is thought to be robust and may be more applicable to other non-highly 
reversible PET ligands with ideal reference tissue. 

Keywords: [^^CjPittsburgh compound B; Alzheimer's disease; Kinetic modeling; PET quantification; Reference tissue; 
Slow dissociation ligand 
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Background 

Various reversible-type radioligands have been developed 
for in vivo neuroreceptor study with positron emission 
tomography (PET). Both arterial blood sampling and long 
dynamic PET scan, up to 120 min, are required for stand- 
ard nonlinear least-squares (NLS) analysis to estimate Ki 
to in the two-tissue compartment four-parameter 
model (4P model): Ki represents the blood-to-brain trans- 
port constant, k2 represents the brain-to-blood transport 
constant, ks represents the first-order association rate con- 
stant for specific binding, and represents the dissoci- 
ation rate constant for specific binding. The k^ represents 
Bmajc'koni where Bmax is maximum receptor density and 
kon is the in vivo association rate constant. Since ks repre- 
sents available receptors for the PET ligand, it is the target 
parameter of major interest in most PET studies. However, 
quantification of ks in the 4P model is often difficult 
because of uncertainty of the /C4 estimate and high correl- 
ation between the ks and /C4 estimates. As surrogate 
parameters for fimax. binding potential and distribution 
volume have been widely used [1-4]. Several reference tis- 
sue methods have also been developed [5-10]. 

Irreversible (enzyme-substrate type) radiotracers ['^^C] 
methylpiperidin-4-yl acetate and propionate have been 
developed for the measurement of cerebral acetylcholine 
esterase activity using PET [11,12]. In this case the two- 
tissue compartment three-parameter {Ki to k^) model 
(3P model) was used to estimate /cs, which is an index of 
acetylcholine esterase activity. In the 3P model, the pre- 
cision of /f3 estimate is usually higher than in the 4P 
model, in spite of shorter PET scan time (40 to 60 min), 
since there is no need of /C4 estimation in the 3P model. 

We have previously defined two mathematical functions, 
the information density function and information func- 
tion, which are useful for model selection and opti- 
mization of scan time in PET [13]. Based on simulations 
using both functions, we proposed a new method (3P + 
method) for quantification of ks for moderately reversible 
ligands. '3P+' means three-parameter model plus short 
PET scan. In this method, the 3P model (/C4 = 0 model) 
was applied to the early-phase PET data (up to 30 to 40 
min) from reversible ligands with moderate k^ (moderately 
reversible ligands). Although the 3P + method was not al- 
ways developed for a specific ligand, the amyloid-binding 
radiotracer [^'^C] Pittsburgh compound B (PIB) was used as 
an example for the moderately reversible ligands (fcj, = 
0.018/min). The 3P + method afforded a more stable ks es- 
timate than the standard 90-min 4P analysis. However, 
there is still the drawback of the necessity for arterial blood 
sampling and radiometabolite analysis, which may restrict 
the widespread use of this method in daily clinical practice. 

In this article, we propose a noninvasive 3P++ analysis 
using [ C]PIB. 3P++ means 3P + analysis plus use of 
a reference tissue for input function. To validate the 



proposed method, the linear correlations of ks estimates 
were evaluated between 40-min 3P++ and 3P + analyses, 
as well as between 3P++ and 90-min 4P analyses in clin- 
ical PET studies. In addition, simulation studies were per- 
formed to explain ks biases observed in the 3P++ analysis. 

Methods 

Theory 

Assumptions in 3P++ analysis 

The following are assumptions used in 3P++ analysis: 

• Assumption 1 (on the nature of radioligand used): We 
apply 3P++ analysis only to moderately reversible or 
nearly irreversible radioligands (ki < G.03/min), but 
exclude highly reversible ligands. [^^C]PIB is an 
example of moderately reversible ligands 

(/C4 = 0.018/min). 

• Assumption 2 (on the duration time of PET scan): 
We use early-phase PET data in the curve fitting. In 
[^^C]PIB, dynamic PET data during 0 to 40 min was 
described well with the 3P model, since the effect of 
the /c4 process on PET data was negligible within 
these early-phase kinetics [13]. 

• Assumption 3 (on the specific binding in the 
reference tissue, /tsJ: Specific binding of radioligand is 
negligible in the reference tissue {ks^ = 0). In [^^C]PIB, 
the gray matter of the cerebellum is usually used as a 
reference tissue for input function [14]. We apply the 
one-tissue compartment two-parameter (Ki, k-i) 
model (2P model) to the reference tissue. 

Working equation for 3P++ analysis 

The working equation for the 3P++ analysis has been 
reported [15]: 



Ctit) = 



m- 



ki + ki 



klr-ko 



k2 + ki 

Rik2{k2 + kj-kir) 
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where Ct{t) is the radioactivity concentration in the tar- 
get tissue and CXt) is that in the reference tissue; /c2r is 
the /c2 in the reference tissue and ® is the convolution 
integral. The rate of tracer penetration into the target 
tissue is obtained as the relative value which is the 
ratio of target Ki to reference K^. 

Clinical PET study 
Human subjects 

Two groups of subjects, a normal control (NC) group 
and an Alzheimer's disease (AD) group, participated in 
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the current study with written informed consent. The 
NC group consisted of 15 healthy subjects (age ranging 
from 48 to 90 years, 66.7 ± 11.5 years (mean ± SD); eight 
males and seven females) without a history of central 
nervous system diseases or psychiatric disorders, and the 
AD group consisted of 15 patients (ages 55 to 85, 68.9 ± 
9.6 years; four males and 11 females) diagnosed as prob- 
able AD according to the criteria of the National Insti- 
tute of Neurological and Communication Disorders, 
Alzheimer's Disease and Related Disorders Association 
[16]. The study was approved by the Institutional Review 
Board of the National Institute of Radiological Sciences. 

Radiochemical synthesis 

[^^C]PIB was synthesized by the reaction of 2-(4'-ami- 
nophenyl)-6-hydroxy-benzothiazole and [ C] methyl tri- 
flate [17]. The product had radiochemical purity greater 
than 95.4%. Specific activity was in the range of 56.3 to 
285.3 GBq/nmol. 

PET scan protocol 

PET images were acquired with a Siemens ECAT 
EXACT HR + scanner (CTI PET systems. Inc., Knoxville, 
TN, USA) with an axial field of view of 155 mm, provid- 
ing 63 contiguous 2.46-mm slices with 5.6-mm tran- 
saxial and 5.4-mm axial resolution. After a 10-min 
transmission scan for tissue attenuation correction, infu- 
sion of ["C]PIB (about 370 MBq in 5 mL for 1 min) 
began. A PET scan in 3D mode was started after the arrival 
of tracer to the brain (approximately 30 s after the begin- 
ning of tracer infusion). The dynamic scans consisted of 19 
frames (3 x 20 s, 3 x 40 s, 1x1 min, 2x3 min, 5x6 min, 
and 5 x 10 min) with the total scan duration of 90 min. All 
data processing and image reconstruction were performed 
using standard Siemens software, which included scatter 
correction, randoms, and dead time correction. 

Region-of-interest delineation 

Region-of-interest (ROI) analysis was performed using the 
PMOD software package (PMOD version 3.2; Technolo- 
gies Ltd., Adliswil, Switzerland). The ["C]PIB PET images 
were co-registered to weighted images in each subject. 
The following 15 ROIs were drawn manually on Ti 
weighted images: frontal, mesial temporal, lateral temporal, 
parietal, occipital, anterior cingulate, and posterior cingu- 
late cortices in both hemispheres as well as the reference 
tissue (gray matter of cerebellum). ROIs were transferred 
to co-registered [^'^C]PIB PET images, and time-activity 
curves (TACs) were obtained in those brain regions. 

Input function measurement 

During PET scan, arterial blood was collected from radial 
artery, starting 6 s (transit delay at the blood sampling site) 
after the beginning of PET scan to 85 min post injection 



(10 x 10 s, 1 x 30 s, 9 x 2 min, 6 x 10 min, and 1x5 min; 
27 samples). Radioactive metabolites were analyzed by a 
radio-thin layer chromatography (TLC) method [12], with 
a TLC-developing solvent (ethyl acetate/«-hexane = 2:1 
vols). The metabolite-corrected radioactivity as well as 
total radioactivity in blood plasma was fitted to a mono- 
exponential saturation ftmction during infusion (0 to 1 
min) and the sum of three-exponential functions after the 
end of infusion (1 to 85 min) [12]. 

4P and 3P + analyses (arterial-plasma input) 

Brain regional TACs were analyzed by the weighted NLS 
method under positive constraint of all /c, with metabolite- 
corrected input function to afford /<i to estimates in 4P 
analysis (scan time of 90 min) and Ki to ks estimates in 
3P + analysis (40 min). Correction was made for blood- 
pool (5%) radioactivity in brain tissue [14]. Custom soft- 
ware operating in IDL software (version 6.0; Jicoux 
Datasystems, Inc., Tokyo, Japan) environment was used 
for the compartment model analysis. 

3P++ analysis (reference tissue input) 

For successful convergence in NLS optimization using 
Equation 1, we fixed k2r to 0.178/min (mean cerebellar /c2 
value by 40-min 3P + analysis; N= 30; SD = 0.034). Based 
on Equation 1 and cerebellar TAC with a fixed ^2r value, 
the time-integral of Cr{t) (the second term on the right 
side of Equation 1) and the convolution integral (the third 
term) were calculated numerically without data inter- 
polation for each scan mid-times during 0 to 40 min, and 
the three parameters Ri, /c2, and k^ were estimated. 

Simulation study 

Generation of error-added TACs for IVIonte Carlo simulation 

The error-free, baseline TACs (19 frames/90 min) simulat- 
ing the target ROI of the NC and AD subjects were gener- 
ated by using the 4P model with parameter set (Ki = 0.180 
mL/g/min, /c2 = 0.180/min, /c3 = 0.018 and 0.036/min for 
the NC and AD subjects, respectively, and = 0.018/min; 
typical values for [^^C]PIB) and averaged (N=20) input 
function of [^'^C]PIB. The reference ROI was the same be- 
tween NC and AD subjects and was generated by using 
the 2P model with parameter set (/<i = 0.180 mL/g/min, 
k2 = 0.180/min) and the same input function as above. The 
error-added TACs for simulation were generated accord- 
ing to the following formula [18]: 

Error-added C, = C, + Rand x ff(Ci), 

where C, is noise-free simulated radioactivity concentra- 
tion at frame number i, Rand is a random number from a 
Gaussian distribution with a mean 0 and variance 1, £ is a 
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scaling factor that determines the noise level, Ati is scan 
duration of frame number /, i, is mid-scan time of frame 
number i, and A is ^^C decay constant. In all Monte Carlo 
simulations, a data set of 100 noise-added TACs was ana- 
lyzed with weighted NLS, using a relative weight w,-: 



w, = constant x 



Ati X exp{-Ati) 



(3) 



Effects of PET noise on 4P, 3P+, and 3P++ analyses 
Five levels of PET noise (0.025, 0.05, 0.1, 0.2, and 0.3; e in 
Equation 2, relative values empirically determined) were 
added to the baseline TACs of the target ROI of the NC 
subjects. From 100 error-added TACs for each PET noise 
level, 100 /c3 values were estimated using 90-min 4P, 40- 
min 3P+, and 3P++ analyses. Coefficient-of-variation (CV) 
of ks was calculated as CV (%) = (SD/mean) x 100. In the 
following simulations, the PET noise was fixed at 0.1. 

Effects of K, change in target ROI on 4P, 3P+, and 3P++ 
analyses 

Simulated target TACs were generated by 4P model with 
five different Ki values (0.12, 0.15, 0.18, 0.21, and 0.24 
mL/g/min) and fixed ^3 (0.018/min) and (0.018/min). 
The value of Ki/k2 was fixed at 1. The range of Ki was 
determined with clinically measured Ki for ['^^CJPIB 
(0.177 ± 0.31 in NC group and 0.168 ± 0.30 in AD group; 
90-min 4P analysis). Reference TAC was the same as 
baseline reference TAC. The bias in 90-min 4P, 40- 
min 3P+, and 3P++ analyses relative to the true k^ 
(0.018/min) was calculated as bias (%) = (estimated /cs/ 
true ^3 - 1) X 100. 

Effects of k2 or fcj change in reference ROI on 3P++ analysis 

In 3P++ analysis, ^3^ was assumed to be 0 and /c2r was 
fixed as an empirical constant. The effects of /c2r or k3^ 
change were investigated as follows. The error-added tar- 
get TACs were generated by 4P model with two different 
^3 values (0.018/min for NC and 0.036/min for AD); other 
parameters were the same as the baseline target TAC. The 
error-added reference TACs were generated by 2P model 
with five different ^2 (0.12, 0.15, 0.18, 0.21, and 0.24/min) 
and fixed /<i values (0.18 mL/g/min). Another set of simu- 
lated reference TACs was generated by 3P model (not 2P 
model) with five different /C3 (0, 0.002, 0.004, 0.006, and 
0.008/min) and fixed Ki (0.18 mL/g/min) and ^2 (0.18/ 
min). The k^ bias in 3P++ analysis was expressed relative 
to 3P + analysis as bias (%) = (3P++ /C3/3P + ^3 - 1) x 100. 

Although /csr was assumed to be 0 in Equation 1, each 
subject may have different k^^ values that deviated from 
0. In simulations to investigate the effect of the individ- 
ual ksr variation on 3P++ analysis, we defined the ks 
value empirically corrected for nonzero k^,- as follows: 



^3 ' = ^3 + fen where ^3 is the k3 estimate of target ROI by 
3P++ analysis and /c3r is the /C3 estimate of reference ROI 
by 3P + analysis (true reference /C3). Bias in 3P++ k^' 
relative to 3P + /C3 was compared with the bias in 3P++ 
ks to 3P + /c3. 

Results 

Goodness of model fits in 3P-I-+ analysis 

Figure lA shows an example of the curve fitting of [''q 
PIB cerebellar TAC data to the 2P model, where a good fit 
is seen during 0 to 40 min after tracer injection. Figure IB 
shows the fits of cerebral cortical TAC data (0 to 40 min) 
to the 3P + and 3P++ models. The goodness-of-fit by 3P++ 
model (reference tissue input) is almost indistinguishable 
from that by 3P + model (arterial-plasma input). Kinetic 
parameters (/<i = 0.161 mL/g/min, /c2 = 0.167/min and 
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Figure 1 Reference and target tissue TACs in [ C]PIB PET. (A) 

Cerebellar (reference tissue) data (open circle) up to 90 min in one 
AD subject and the fit of 40-min data to the 2P model (solid line). 
(B) Cerebral cortical (target tissue) data (open circle) in the same 
subject and the fits of 40-min data to the BP model with 
arterial-plasma input (3P 4- analysis; solid line) or reference tissue input 
(3P-H- analysis; dashed line). The dotted lines in (A) and (B) indicate 
the extension of the solid line from 40 to 90 min. 
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/C3 = 0.015/min) were estimated in 3P + analysis and Ri = 
0.897, /C2 = 0.158/min and /cs = 0.011/min in 3P++ analysis. 

Intra-subject correlation 

Figure 2A is an example of the intra-subject correl- 
ation between 40-min 3P + (X-axis) and 3P++ (F-axis) 
analyses, where the /cs values of 15 ROIs, including the 
cerebellum (reference tissue in 3P++ analysis) from one 
particular NC subject or one particular AD subject, are 
shown. The regression lines and the coefficients of de- 
termination are r= 0.845X- 0.006 {/ = 0.972) for the 
NC subject and r= 0.655X- 0.004 (r^ = 0.982) for the 
AD subject. Cerebellar /cs values for both subjects are 
naturally calculated to be 0 in the 3P++ analysis. The 
slopes of the regression lines indicate the presence of 
negative bias in the 3P++ against the 3P + analysis. 

Figure 2B shows the /cs correlation between 90-min 4P 
(X-axis) and 40-min 3P++ (^-axis) analyses in the same 
subjects. The regression lines are r=0.590X- 0.005 {/■ = 
0.953) for the NC subject and r= 0.338X+ 0.000 (r^ = 
0.907) for the AD subject. When the cerebellar data 
(X= 0.008, >'= 0.000) was removed from calculation for 
the AD subject, the regression line became y=0.295X- 
0.002 with slightly larger (0.935; not shown in the 
figure). The slopes of the regression lines show that /cs bias 
in 3P++ against 4P analysis is larger than that against 
3P + analysis. 

Inter-subject correlation 

Figure 3A shows an example of the inter-subject /cs cor- 
relation, where /cs values for the left lateral temporal 
cortex from 30 subjects (15 NC +15 AD) are compared 
between 40-min 3P + (X-axis) and 3P++ (7- axis) ana- 
lyses. The regression lines are r=0.461X- 0.001 (r^ = 
0.739) for all 30 subjects, Y= 0.178X+ 0.000 (r^ = 0.151) 



for the NC group alone, and r= 0.286X+ 0.003 (r^ = 
0.411) for the AD group alone; the latter two lines are not 
shown in the figure. The slopes of the regression lines also 
indicate the presence of negative biases in 3P++ against 
3P + analysis. 

Figure 3B shows the inter-subject correlation of left 
lateral temporal between 90-min 4P (X-axis) and 40- 
min 3P++ (y-axis) analyses, where the regression line is 
r=0.225X+ 0.000 (r^ = 0.711) for all subjects. The lines 
of Y= 0.090X+ 0.001 (r^ = 0.122) for the NC group alone 
and y=0.135X+ 0.005 (r^ = 0.513) for the AD group 
alone were also calculated. The slopes of the regression 
lines show larger negative k^ biases in 3P++ against 4P 
analysis than that shown in Figure 3A. The results in 
other cerebral regions were essentially the same as those 
in the left lateral temporal cortex. 

Simulation on the effects of PET noise on ks CV 

Figure 4 compares the noise sensitivity of k^ estimates 
among the 90-min 4P, 40-min 3P+, and 3P++ analyses. In 
all three analyses, the k^ CVs increased as the PET error 
became larger. The k^ CV in 3P++ analysis was compar- 
able to that in 3P + analysis and lower than that in 4P ana- 
lysis; for example, k^ CVs at 0.1 of noise level were 6.6% in 
3P++, 7.0% in 3P+, and 11.4% in 4P analyses. 

Simulation on the effects of target change on kj bias 

Figure 5 shows the effects of Ki change in the target 
ROI on the /cs biases in the 90-min 4P, 40-min 3P+, and 
3P++ analyses. The 4P analysis remained almost bias- 
free (+0.6%) within /Ci from 0.12 to 0.24 mL/g/min. 3P 
+ and 3P++ analyses showed larger negative biases 
(-33% to -34% bias in 3P + and -33% to -35% bias in 
3P++) compared with 4P analysis. Although 3P++ ana- 
lysis showed slightly larger /fs bias than 3P + analysis 
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Figure 3 Inter-subject correlations of left lateral temporal in ["CjPIB PET. The results in 40-min 3P++ (/-axis) vs. 3P + (X-axis) analyses 
(A) and 40-min 3P4-+ vs. 90-min 4P analyses (B) with 30 subjects (IS NC, open circle; 15 AD, closed circle) are shown. 



when Ki was low (0.12 mL/g/min), bias in 3P-i--t- ana- 
lysis was almost the same as 3P -t- analysis. 

Simulation on the effects of change on 3P-i~h analysis 

In 3P-i--t- analysis (Equation 1), /c2r was fixed at 0.178/ 
min, though k2r was not always the same among subjects 
(CV = 19%). Figure 6 shows the effects of individual ^2r 
change in 40-min 3P-H- analysis. When /c2r was equal to 
the fixed value (0.18/min), 3P++ analysis was bias-free, 
relative to 3P -1- analysis. However, when /c2r was different 
from the fixed value, 3P-t--h analysis showed a negative ks 
bias relative to 3P + k^. The /c2r effects were similar 



between NC ROI (/fj = 0.018/min) and AD ROI (/C3 = 
0.036/min); for example, the biases were -14.1% for NC 
and -12.1% for AD at /c2r = 0.12/min and -14.1% for NC 
and -11.3% for AD at /c2r = 0.24/min. 

Simulation on the effects of V.^, change on 3P-n- analysis 

In 3P+-H analysis we assume that k^^ = 0, that is, specific 
binding is negligible in the reference tissue. However, in 
all subjects examined, this assumption did not hold: the 
kzr values in 40-min 3P + analysis were 0.008 ± 0.004/min 
in the AD group, 0.007 ± 0.002/min in the NC group, and 
0.007 ± 0.003/min in the AD + NC group. 
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Figure 4 Effects of PET noise on CV of k^. The results in 40-min 
3P++ (open triangle), 40-min 3P + (open square), and 90-min 4P 
(open circle) analyses are shown. Five different PET noises (0.025 to 
0.3) were added to the ["C]PIB baseline TACs of the target ROI of 
the NC subjects. CV of was calculated from 100 kj, estimates as 
CV (%) = (SD/mean) x 1 00. 
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Figure 5 Effects of K-t change in tlie target region on k^ bias. 

The results in 40-min 3P++ (open triangle), 40-min 3P + 
(open square) and 90-min 4P (open circle) analyses are shown. 
Simulated target TACs were generated by 4P model with five 
different values (0.12 to 0.24 mL/g/min). The k^ bias was 
calculated as bias (%) = (estimated /ts/true /<3 - 1)x 100. 
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0.12 0.15 0.18 0.21 0.24 

reference k2 (/min) 

Figure 6 Effects of k2 change in the reference region on bias 
in 40-min 3P++ analysis. Simulated target TACs were generated by 
4P model with two different values {0.018/min for NC, open circle; 
0.036/min for AD, closed circle). Simulated reference TACs were 
generated by 2P model with five different k2 values (0.12 to 0.24/ 
min). The bias in 3P++ analysis was expressed relative to 3P + 
analysis as bias (%) = (3P++ k-^/iP + - l)x 100. 



Figure 7 shows the effects of individual /cs^ change (0 to 
0.008/niin) on 40-min 3P++ analysis. When ks, was 0, 3P++ 
analysis was bias-free, relative to 3P + analysis. The 
biases (negative biases) increased as /csr increased: -38% 
for NC and -27% for AD at ks^ = 0.004/min and -70% for 
NC and -48% for AD at ks^ = 0.008/min. The NC ROI 
= 0.018/min) showed larger biases than the AD ROI 
{ks = 0.036/min). Figure 7 also shows the results of the 
simulation study on the relationship between 3P++ /cs' 
and 3P + ks, where 3P++ /cs was empirically corrected with 
individual ks^- In this case, negative bias in 3P++ ks' was 
significantly decreased compared to that in 3P++ /cs; for 
example, bias was decreased from -70% to -7% for NC, 
and from -48% to -15% for AD at ks, = 0.008/min. 

Figure 8 shows the correlation between 3P++ ks' and 
BP + k^ using the same data as in Figure 3A, where 3P++ 
/c3 in Figure 3A was replaced by 3P++ ks'. The regres- 
sion line was r= 0.678X+ 0.003 (r^ = 0.975) for all sub- 
jects, where X=3P + k^ and Y= 3P++ k^'. The lines of Y 
= 0.798X+ 0.002 (r^ = 0.897) for the NC group alone and 
r=0.620X+ 0.004 (r^ = 0.960) for the AD group alone 
were also calculated. The determination coefficient was 
increased by this correction from 0.739 to 0.975. The 
slope of the regression line was also increased from 
0.461 (Figure 3A) to 0.678 (Figure 8), which showed the 
reduction of negative bias in 3P++ analysis. 

Discussion 

Theoretical basis and merits of 3P++ analysis 

The previous 3P + analysis allowed for estimating /cs of 
moderately reversible ligands, where the 3P model was 



+ 

CL 
CO 



+ 
+ 
Q. 

CO 



-20 



-40 



? -60 



^^T^^^A^-^ . . . NC /cV 


\\ ^ — ^ 




AD/C3' 




AD/C3 


1 1 1 


NC /C3 



-80 

0.000 0.002 0.004 0.006 0.008 0.010 

reference (/min) 

Figure 7 Effects of change in the reference region on ^3 bias 
in 40-min 3P++ analysis. Simulated target TACs were generated by 
4P model with two different k-^ values (0.018/min for NC, open circle; 
0.036/min for AD, closed circle). Simulated reference TACs were 
generated by 3P model with five different k^ values (0 to 0.008/min). 
The /C3 bias in 3P++ analysis was expressed relative to 3P + analysis 
as bias (%) = (3P++ k^/iP + /cj - 1) x 100. Effects on bias in 3P++ k^' 
relative to 3P + k^ are also shown (NC, open triangle; AD, closed 
triangle), where 3P++ k-^' was calculated as (3P++ k^') = (3P++ 
/C3) + reference k-^. 



applied to early-phase (up to 30 to 40 min) PET data 
with arterial input function [13]. It was reported that 
when the 3P model was applied to 60-min PET scan data 
from ['X]PIB (^4 = 0.018/min) as a moderately revers- 
ible ligand, only a poor model fit was obtained [19]. Pre- 
vious simulation studies on [ C]PIB using information 
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Figure 8 Inter-subject correlation of left lateral temporal k^ in 
["C]PIB PET. The result in 40-min 3P++ (/-axis) vs. 3P + (X-axis) 
analyses with 30 subjects (15 NC, open circle; 15 AD, closed circle) is 
shown. The k-^ estimates were empirically corrected as (3P-H- k{) = 
(3P-H- y + (individual cerebellar k^ by 3P+). 
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density theory suggested that scan time reduction to 40 
min would be necessary to obtain a good fit to the 3P 
model [13]. 

When 3P + or 3P++ analysis can be applied to a ligand, 
such ligand is specified as a moderately reversible ligand. 
This applicability is determined by the information func- 
tion curves of 1(3 and 1(4, [13], and thus is dependent on 
the scan time as well as ks and 1(4 values of the ligand in 
a ROI. Differentiation of a moderately reversible ligand 
from general reversible ligands is somewhat arbitrary, 
though we conveniently defined this with the 1(4 value 
(<0.03/min) in this study. 

In the present study, the 3P + plasma input model was 
extended to the 3P++ reference tissue input model. The 
3P++ analysis has three merits over previous methods. 
First, the PET scan time is short, usually less than 40 min, 
which may be important in PET studies with elderly or de- 
mented subjects. Secondly, the target parameter can be 
isolated from the other model parameters. Thirdly, neither 
arterial cannulation nor labor-intensive measurements of 
labeled metabolites are required. 

One of the conventional models for the estimation of 
binding of ["^^CIPIB is the Logan plot analysis [2], which 
employs data of long duration (more than 60 min). Non- 
invasive Logan analysis (distribution volume ratio) [6] re- 
quires late-phase (equilibrium-phase) PET data, whereas 
late-phase data are not necessary for 3P++ analysis. In the 
noninvasive Logan model or simplified reference tissue 
model [8], the Ki-to-k2 ratio in the target and reference tis- 
sues is assumed to be equal. 3P++ analysis does not require 
such an assumption. Since 3P++ analysis is a kind of 
irreversible-model analysis, Ki (Ri) and ks can be independ- 
ently estimated (^2 must be frxed to a certain constant). 

Noise sensitivity of 3P++ analysis 

Loss of PET data in short-scan 3P++ and 3P + analyses 
might be considered to deteriorate the precision of the 
/c3 estimate. In the present simulation for noise sensitiv- 
ity, ks CV values in 40-min 3P++ and 3P + analyses were 
lower than (almost three fifths of) that in 90-min 4P 
analysis (Figure 4), which was in accordance with the 
previous report [13]. It is considered that the loss of 
PET data may be compensated for by the reduction in 
the number of free parameters from four in the 4P 
model to three in the 3P + and 3P++ models. 

Ki effect on 3P++ analysis 

In the Ki simulation, the stableness of k^ estimation in 
changes of cerebral blood flow was investigated. The mag- 
nitudes of ks bias were independent of the Ki change, ran- 
ging from 0.12 to 0.24 mL/g/min, in 3P++, 3P+, and 4P 
analyses (Figure 5). The 3P++ as well as 3P + and 4P ana- 
lyses were less affected by Ki, which is owing to the cap- 
ability of isolating the k^ estimation. The 40-min 3P + 



analysis showed -33% /cs bias relative to 90-min 4P ana- 
lysis, which is in accordance with the previous report [13]. 
In this Ki simulation, 3P++ k^ showed negligible bias rela- 
tive to 3P + ki. These results suggested that in 3P++ ana- 
lysis, the effects of ignoring vascular volume as well as 
numerical integration error due to discrete time points 
were not significant. 

Causes of negative bias in 3P+-I- analysis 

Firstly, the /cs bias in 3P++ analysis originates from 3P 
model approximation. Our previous simulation study 
[13] showed that the 3P + analysis with 28-min scan had 
large negative ki bias relative to 4P analysis with 90-min 
scan; for example, there was about -22% to -24% bias to 
true /c3 (4P /C3) ranging from 0.01 to 0.04/min including 
NC and AD /C3. 3P++ analysis showed further negative 
ks bias relative to 3P + analysis due to the following two 
reasons. 

Secondly, the bias is due to individual k2r change from 
the fixed value in Equation 1. In 3P++ analysis, we also 
assumed that k2 in the reference tissue was constant and 
was fixed at 0.178/min, which was the average /c2 value 
with the 3P + model. In simulation, negative ks bias was 
predicted when /c2r was larger or smaller than fixed k2 
(Figure 6). Each subject in the NC and AD groups had 
different /c2 values in the reference tissue, and it is con- 
sidered that such biological variance as for reference tis- 
sue may result in a negative ks bias in 3P++ analysis, 
relative to 3P + analysis for [^'^C]PIB. 

Thirdly, the bias is due to the discrepancy between the 
model assumption and the actual reference ROI. The basic 
assumption (assumption 3) in 3P++ analysis is /csr = 0. The 
working equation of 3P++ analysis (Equation 1) is derived 
under this assumption, and reference ki is naturally calcu- 
lated to be 0. However, in 3P + analysis with [^^C]PIB, the 
cerebellum showed nonzero ^3 (0.007 ± 0.003/min in all 30 
subjects). Thus, 3P++ ks is expected to be underestimated. 
Simulation studies showed that 3P++ analysis was bias-free 
for ideal reference with zero k^ and that k^ bias became lar- 
ger as kir increased (Figure 7). When ki was replaced by 
ki', negative bias was significantly decreased in the simula- 
tion (Figure 7), as well as the slope of the regression line 
between 3P++ and 3P + analyses being increased from 
0.461 (Figure 3A) to 0.678 (Figure 8), which also suggested 
that nonzero ki^ caused underestimation of 3P++ ki. 

Correlation of ks between 3P++ and 3P + analyses 

Strong intra-subject ki correlation was shown between 
3P++ and 3P + analyses, and the rank-order of ki was al- 
most the same between the two analyses (Figure 2A), 
suggesting the stability of both 3P++ and 3P + analyses. 

The inter-subject ki correlation (r^; Figure 3A) was 
significantly lower than the intra-subject correlation 
(Figure 2A). Such a lower inter-subject ki correlation 
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can be partly explained by the sample variance of cere- 
bellar ks. In order to explain this, ks' was calculated for 
each subject. When was replaced by k^', the deter- 
mination coefficient between 3P++ and 3P + analyses 
was increased from 0.739 (Figure 3A) to 0.975 (Figure 8); 
the latter is comparable to of the intra-subject Ats cor- 
relation (0.982; Figure 2A). 

Such an estimation of parameter /cs ' is not always prac- 
tical, as 3P + analysis with arterial input function is neces- 
sary for individual cerebellar ks estimation. However, these 
results suggest that the lower in the inter-subject correl- 
ation compared with the intra-subject correlation is due to 
the sample variance of cerebellar /cs and that 3P++ analysis 
itself is robust, as far as the reference is ideal. 

Practically, the use of mean k^,- may be meaningful. 
When target ks is empirically corrected as corrected /cs = 
estimated ka + mean cerebellar k^, the absolute bias in tar- 
get /c3 would decrease. However, the precision of target k^ 
would not necessarily be improved owing to the variance 
of individual /csr- 

In addition to the nonzero effect of k^p inter-subject 
variation of /c2r from the fixed value (/c2 = 0.178/min) 
may also produce individually different ks bias in 3P++ 
analysis, resulting in lower inter-subject k^ correlation 
between 3P + and 3P++ analyses. 

Limitations of 3P-I-+ analysis 

When 3P++ analysis was applied to ["qPIB as an ex- 
ample of moderately reversible ligands, a somewhat lower 
inter-subject ks correlation (r^ = 0.739 or 0.711; Figure 3A 
or Figure 3B) was shown between the 3P++ and 3P + or 
4P analyses, respectively, across a /cs range including NC 
and AD (3P + /cg, 0.004 to 0.040/min). The rank order of 
3P++ A3 also differed considerably from 3P + /cs or 4P /C3. 
These results were mainly due to nonzero ^31 and the 
sample variance of both /c2r and /cs^ as described above. 
The negative ks bias (3P++ vs. 3P+) was larger in NC ROI 
(-70%) than in AD ROI (-48%) when A:3r = 0.008/min 
(Figure 7). The previous report showed that the difference 
in /c3 bias (28-min 3P + vs. 90-min 4P) was small between 
NC ROI (-23%) and AD ROI (-24%) [13]. Therefore, the 
/c3 value in 3P++ analysis may be somewhat underesti- 
mated in the ROI with lower amyloid deposition com- 
pared to 3P + or 4P analysis. 

In [^^'^CJPIB PET, 3P++ analysis may be inadequate 
for inter-subject A3 comparison and useful only for 
intra-subject (inter-ROI) comparison or pre- vs. post- 
comparison in the same subject. 3P++ analysis would 
be more suitable for such reversible ligands that have 
moderate A4 and reference tissue without specific binding. 

Conclusions 

The 3P++ analysis is a A3 estimation method for mod- 
erately reversible PET ligands with a short scan time 



such as 40 min and without arterial blood sampling. 
Although the applicability of 3P++ method to ["C]PIB 
PET may be restricted to intra-subject comparison, 
3P++ analysis itself is robust. The 3P++ method would 
be useful for PET study with non-highly reversible li- 
gands, as far as the reference tissue without specific 
binding is available. 
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